The time evolution of coherent atomic system and probe light in an EIT medium 
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The adiabatic solutions of Maxwell-Bloch equation gov- 
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evolution of the density matrix elements of the EIT system 
and the probe light is thus investigated by using the adiabatic 
approximation formulation and the slowly varying envelope 
condition. 
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I. INTRODUCTION 

Recently, many theoretical and experimental investi- 
gations show that the control of phase coherence in mul- 
tilevel atomic ensembles will give rise to many novel and 
striking quantum optical phenomena, such as the coher- 
ent population trapping [1], laser without inversion and 
electromagnetically induced transparency (EIT) [2-4] , in 
the wave propagation of near-resonant light. According 
to the theoretical analysis of multilevel atomic phase co- 
herence, the requirement of the occurrence of EIT is such 
that the strength of coupling light is much stronger than 
that of probe light [4,5]. Under this condition, the EIT 
atomic vapor allows the probe light to propagate with- 
out dissipation through the medium. Due to its unusual 
quantum coherent character, the discovery of EIT has so 
far led to many new peculiar effects and phenomena [4], 
some of which are believed to be useful for the devel- 
opment of new techniques in quantum optics. More re- 
cently, the physical effects associated with EIT observed 
experimentally include the ultraslow light pulse propa- 
gation [6,7] and light storage [8,9] in atomic vapor, and 
atomic ground state cooling [10]. 

With the development of the quantum information 
(quantum computation, quantum communication and 
quantum measurement), the search for new ways to ma- 
nipulate photon states becomes increasingly important 
[11,12]. In the area of EIT, historically, such a possibil- 
ity might first be considered by Ling et ai, who inves- 
tigated the "electromagnetically induced grating" in ho- 
mogeneously broadened media [13]. In this work there is 
a strong coupling standing wave, interacting with three- 
level Lambda-type (or cascade-type) atoms. This can 



diffract a weak probe field, which propagates along a 
direction normal to the standing wave, into high-order 
diffractions. As stated by Arve et al, many researches 
of EIT including the freezing light pulses [14,15], time- 
dependent group velocity [16,17], delayed probe light [8], 
stopping a pulse and restarting it in the opposite di- 
rection (creating a time reversed version of the origi- 
nal probe pulse) [18] have been limited to electromag- 
netic propagation in one dimension only [19], namely, 
the investigators of these studies did not consider how 
to change the direction of the probe light. For this rea- 
son, Arve et al. lifted the restriction to the one dimen- 
sion and considered the Maxwell-Bloch dynamics in two 
or more dimensions of three-level EIT system, by using 
the adiabatic approximation and the slowly varying en- 
velope approximation. In their treatment, however, they 
did not take into account the decay terms in the dy- 
namical equations. Since in a typical EIT experiment 
[20,21], the spontaneous decay rate of the excited state 
of the three-level systems may be of the same order of 
magnitude as that of the Rabi frequencies of the probe 
and coupling optical fields, the contribution of the de- 
cay terms in the Maxwell-Bloch dynamics should also be 
taken into consideration. In the present paper, we obtain 
the adiabatic solutions of the Maxwell-Bloch equations, 
i.e., the expressions for the time evolution of both the 
probe light and the off-diagonal density matrix elements 
of the atomic system. The formulation presented here 
can apply to the study of the variation of the propaga- 
tion direction (diffraction) of the optical fields in the EIT 
medium. 

II. MAXWELL-BLOCH EQUATIONS 
GOVERNING THE THREE-LEVEL EIT SYSTEM 

Consider a three-level A-type atomic ensemble inter- 
acting with two resonant laser beams, the Rabi frequen- 
cies of which are denoted by Q c and Q p , respectively. In 
such a A-type atomic system, levels |c) and \b) are the 
ground states, and a) the excited state. The laser beam 
which couples levels \a) and \b) is called the probe light 
(fip). Another laser beam is termed the coupling light 
(f2 c ), which couples levels \a) and \c). The schematic di- 
agram for the above A-type system is depicted in Fig. 
1. 
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If we apply the adiabatic approximation (i.e., 
\dp/dt\ <C |fi c |) to the three- level system, then the Bloch 
equation can be rewritten in the form 



= Im (0*/? ah + n*p ac ) - laaPaa, 







[flcPcb 



Op {Pbb - Paa)] - "fabPab, 
= f l^pPbc + O c (Pec - Paa)} ~ JacPac, 
p bb = Im(O p p ha ) - "IbbPbb, 
Pbc = 5 (0*/9 ac - QcPba) - IbcPbc, 
= lm(n c p ca ) -JccPcc, 



(1) 



where dot denotes the time derivative. In order to let 
level b) into a dark state (trapped state) and thus realize 
the transparency effect for the probe light, the magnitude 
of Rabi frequency |f2 c | should be much greater than |O p |. 
Generally speaking, 0*p a0 in the first equation in (1) 
can be ignored compared with the term 0*p ac . It follows 
from the first and last equations in (1) that the relation 
laaPaa + IccPcc = can be satisfied, and consequently 
the ratio p aa / ' p cc = -j cc /-f aa is obtained. 

Assume that the initial state of the excited level \a) 
of the EIT medium under consideration is nearly empty, 
and the two ground states |c) and \b) are occupied ac- 
cording to the dark state conditions, i.e., pbb — 1 and 
p aa — 0. This, therefore, means that the term pbb — Paa 
in the second equation in (1) approximately equal unity. 
So, we can obtain the following relation between the den- 
sity matrix elements p c b and p a b 



Pcb 



Op + 2i-f a bPab 
fir. 



(2) 



In the meanwhile, it follows from the fifth equation in (1) 
that the relation between pb a and pb c is 



Pba 



2 * jmPbc 



IbcPbc 



(3) 



where pb a and pb c are the complex conjugations to p a b 
and p c b, respectively. Note that Eqs. (2) and (3) consti- 
tute a set of coupling equations. In order to see how pb c 
and pba evolve in the presence of the probe and coupling 
lasers, we should first obtain their respective equations. 
With the help of Eqs. (2) and (3), one can arrive at 



d 

dl Pbc + 



and 



dt Pba 
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(4) 



27a6 



(5) 



As in the adiabatic approximation, the Rabi frequency 
O c of the coupling laser can be considered a constant in 
time, one can derive pb c (t) from Eq. (4), i.e., 



Pbc(t) 



-At 



pV ' e xt dt' + p hc (0) 



(6) 



with A = 7&c + fl*n c /(4:-fab)- In the meanwhile, one can 
obtain the expression for pb a {t) from Eq. (5) 



Pba(t) = e 



— e -At 
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+Pba(0)e 



-At 



(7) 



It should be noted that although the Rabi frequency fl c 
is a time-dependent quantity, in the adiabatic case, it can 
be viewed as a constant number. Such a result will be 
theoretically validated by using the Maxwellian equations 
in the following. Thus, the parameter A in Eqs. (6) and 
(7) is considered a constant number. 

In the above, we consider the Bloch equation governing 
the three-level EIT atomic media. In what follows, we 
will discuss the Maxwellian equations 




lc i(k c 
2 o c e 



-OJ c t) 



e tWVac-, 
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2 6 P e e dt?rab 



(8) 



of the coupling and probe lasers in an EIT medium, 
where N and eo denote the atomic number den- 
sity (total number of atoms per volume) and the 
electric permittivity in vacuum, respectively. Note 
that here p a c{t) and p a b(t) can be rewritten as fol- 
lows p ac (t) = p ac (t)exp [\ (ujact - k ac • r)] , p ab (t) = 
Pabit) exp [\ (uj ab t - k ab ■ r)] . By using the slowly vary- 
ing envelope approximation (i.e., the second order deriva- 
tives of the envelopes £ CjP are negligible compared with 
the other terms on the left-handed sides of the above 
Maxwellian equations), Eq. (8) can be reduced to the 
first order differential equation 



dt° c 



1 Pac? 
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where 
Pab(t) 



= e 



f — 

J a -2ijab 



dt 



z xt 'dt' 



+Pab(0)t 



-At 



(10) 



Here p ca and pb a stand for the transition dipole matrix 
moments. Note that because of the initial condition (e.g., 
the excited state \a) is nearly empty and, moreover, the 
population probability of level |c) is negligibly small), 
at least for the case of investigating the transient optical 
properties of EIT media, p ac can be taken to be zero [19]. 
Thus the wave equation of the coupling light £ c takes the 
form of a source-free field equation. Since such an equa- 
tion is simple, we will not further consider it. Instead, in 
the next section we will concentrate our attention on the 
wave equation of the probe light 
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Ck p • VS2p + — Sip = I ^ p ah , 



(11) 



where the coefficients involved are defined as 



where fl p — £ p -p a b/h with p a b being the transition 
dipole matrix moment of ab transition. 

It is apparently seen that if we have solved the solu- 
tion il p of Eq. (11), then according to (6) and (7), the 
information on the time evolution of the medium polar- 
ization (pbc and pba) can therefore be easily obtained. 
Such a study can apply to the investigation of the adia- 
batic storing of a light pulse by the EIT mechanism [19]. 



III. TIME EVOLUTION OF THE PROBE LIGHT 

The solution f2 p of Eq. (11) takes the form 

{l p (r,t)=f(r-%j\ e (t')dt'y (12) 

where k p is a unit vector defined as k p = k p /|k p |, and 
v g (t) the group velocity of the probe light in the EIT 
medium. Here / is a certain function. It is readily veri- 
fied that the spatial derivative of fl p can be expressed in 
terms of the time derivative of f2 p , i.e., 



k p .VQ p = -l|0 p . 



Thus, Eq. (11) can be rewritten as 



1 I — = l— r 

e h 



Vc, ) dt 



Pab- 



(13) 



(14) 



Further calculation (i.e., calculating the time derivative 
of the above equation) yields 

Vg 9 ( C \ d 2 UpNpbaPab d 

^d-t^ + { 1 --J^ = l ^^dt Pab - (15) 
In accordance with Eq. (5), we have 



d_ 

dt 



—i>ab = ~ x P ab + ^~ b ('^y-p • " ; ' - ! P ) ■ : > ! >' 



Substitution of Eq. (14) into (16) yields 

d ( .UJpNpbaPab^ 1 ( c\ d 
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+ 7. ( — + 7fcc^p I ■ 



27a6 \dt 



(17) 



It follows from Eqs. (15) and (17) that the Rabi fre- 
quency of the probe field, fi p , agrees with the following 
equation 



q2 q 



(18) 



r - 1 v s 



2-y ab e h(l-±) 



UpNpbaPab 



(19) 



In this paper, we consider the homogeneous EIT me- 
dia only, where the group velocity of the probe light 
is independent of time, i.e., v g = 0. Thus the co- 
efficients £ and cr in Eq. (18) are constant. Such a 
choice will simplify the problem under consideration, 
since the expression (12) can be reduced to the simple 
form f2 p (r, t) = f (r — k p u g ^ . 

Apparently, the general solution of Eq. (18) may be of 
the form 

ft P (r,t) = Q p+ (r)exp(?7 + i) + Q p _(r)exp(?7_i), (20) 

where f2 p ±(r) is defined 

O p± (r) = fi p± (0,0)cxp(crk p -r) (21) 

with fi p ±(0, 0) = fi p ±(r = 0, t = 0). Here a is a parame- 
ter that characterizes the field amplitude shape. Keeping 
k p • k p = 1 in mind, we rewrite the time-dependent factor 
exp (r)±t) in Eq. (20) as follows 



exp (rj±t) = exp 
Hence Eq. (20) is rewritten 
fi p (r, t) = O p+ (0, 0) exp |crk p 



-ak p .k p (-^)i 



(22) 



+fi p _ (0, 0) exp {ak p • [r - k p (- t] } . (23) 

Comparing the expression (23) with (12), one can arrive 
at the group velocity of the probe light 



a 



(24) 



where r\± = (— ( ± \/ C 2 — 4<r)/2, which is expressed in 
terms of the parameters £ and <r of Eq. (18). Note that 
here v g is constant, according to (19), we have 



C = A + 



VpNpbaPab 



(25) 



Further calculation shows that the group velocity v g 
satisfies the following quadratic equation 



a 2 v 2 - av g ( + <r = 0. 

For convenience, we set ( = A + [3(1 — c/v g ) 
/?7bc(l - c/-Ug)~\ where 

^ _ U p NpbaPab 
2jabe Tl 



(26) 



(27) 
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and rewrite Eq. (26) as 

2 13 + A + ere Acre + /?7 fcc 



= 0, 



the two roots of which take the form 



(28) 



v g ± 



/3 + X + ac± y/(/3 + A + acf - 4 (Acre + /3 7bc ) 



2ct 



(29) 



Note that in a typical EIT experiment, the relationships 
between the parameters -jab, 76c and £! c are as follows: 
lab ~ ft c ~ 10 8 a" 1 , 76c 0.0l7 Qb and w p ~ 10 15 
s _1 [20,21]. Under these conditions, one can verify that 
(P + A + ere) 2 3> 4 (Acre + ft'jbc)- Thus, the group veloc- 
ity, w g _ , of the probe laser in the EIT medium is reduced 
to 



{0 + A + ere) 



v„- = 



4(Acre+/3 7be ) 



2(7 



Acre + /3 7bc 
a (J3 + A + crc) 



£ ?i -|-" c li, 



-C. 



(30) 



It is worth pointing out that the above result (30) is self- 
consistent, as it can be validated by using the formula of 
group velocity, v g = c/(n + ujdn/dto), where the optical 
refractive index n(uj) = y/Y+xW)- Here, the electric 
susceptibility for the probe laser with a mode frequency 
cj in a steady EIT system is of the form [22] 



xM = - 



Np b aPab 



LO - LO ab + ij bc 



£0?i (u/ - U) ab + lj bc ) (u> - UJab + ilab) 



(31) 



In addition, it should be noted that the group velocity 
v g+ (v g+ = —rj- /a) should be ruled out, since it does 
not satisfy the relation ry_ = (— £ — ( 2 — 4<r)/2. So, 
the only retained group velocity of the probe laser in the 
EIT medium is w g _. 

Eqs. (12) and (13) shows that the adiabatic solutions 
obtained in this paper can deals with the problems of 
two- or three-dimensional wave propagation, including 
diffraction. As is shown in Eq. (30), the coherent control 
of the probe light can be realized: specifically, the change 
of fl c will has an influence on the probe group velocity 
(and hence the parameters £ and c). Thus, according to 
Eqs. (12), (13) and (18), the variation of the direction of 
propagation of the probe light may arise in response to 
the change of the field strength of the coupling laser. 



the time evolution of the off-diagonal density matrix el- 
ements of the EIT atomic system. By inserting the ex- 
pression (20) into Eqs. (6) and (7), one can obtain the 
explicit expressions for p bc and pb a in the time-evolution 
process, i.e., 

Pbc(t) = -^Mfa+ + A)- 1 ^^' - e" A *) 

^lab 

+ ( V - + A)- 1 ^_(e"-* - e- xt )} + Pbc (0)e- Xt , 
Pba(t) = + lbc)(v+ + A)- 1 ^ + (e" + t - e- At ) 



+ (V-+lbc)(V-+^r 1 ^ P -(e n - t -e 
+pba(0)e' 



-Xt\ 



-At 



(32) 



In a typical EIT experiment (N = 10 18 m 3 [21]), one 
can have £ ;§> So, r) + , i.e., (— ( + \/ ( 2 — 4?)/2 will be 
reduced to — <r/C that is less than jbc- In the meanwhile, 
A = 76c + ^^c/(47ab), which is two or three orders of 
magnitude larger than 7 bc in a typical EIT experiment 
[20,21]. This, therefore, means that when the time t is 
taken the several Rabi oscillation periods (~ 10~ 8 s), 
e »?+t ~ i an d e -At ^ q Thus, p ba m (32) approximately 
equals 



Pba 



1 



^ilab A 



p+ 



or 



Pba 



76c 



2 7a6 7 6c + 



O* 



(33) 



(34) 



According to Eq. (20), we have Q p (r,t) ~ f2 p+ (r), 
O p (r, t)* ~ 0* + (r). Thus, by using the definition of the 
dipolc-transition electric susceptibility 



X 



2N\ Pab ^ 



-Pab, 



(35) 



one can arrive at the electric susceptibility at probe fre 
quency, i.e., 



xK>) 



^|Pa6 



«76c 



e 0^ "fablbc + 



(36) 



It should be noted that this expression is consistent with 
(31), so long as the relation ui = u) ab (the resonance con- 
dition) is inserted into (31). It is thus believed that 
the obtained adiabatic solutions (32) of Maxwcll-Bloch 
equation governing the three-level EIT medium may be 
self-consistent. This set of solutions can be utilized to 
consider the information storage in a atomic vapor. 



IV. OFF-DIAGONAL DENSITY MATRIX 
ELEMENTS 

In the previous section, we considered the wave propa- 
gation of the probe laser. In this section, we will discuss 



V. MAGNETIC-DIPOLE TRANSITION AND 
PERMEABILITY FOR THE PROBE LIGHT 

In the previous sections, when treating the propagation 
of the optical fields in the three-level system, we did not 
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consider the magnetic-dipole transition, but the electric- 
dipolc transition. It is believed that the contribution of 
the magnetic-dipole transition to the wave propagation 
of the probe light should be taken into account in the 
case of a stronger coupling field in the EIT medium. Ac- 
cording to the appendix to this paper, the density matrix 
elements p a b and p c b satisfy the following matrix equation 



d_ 

dt 



Pab 
Pcb 



Tafc 

in? 



be 



Pab 
Pcb 



2 1S P 



(37) 



Here, r ab = -f ab + iA ab and T hc = j bc + «(A a6 - A ac ), 
where 7^ and jbc represent the spontaneous decay rate 
of level | a) and the dephasing rate (nonradiative decay 
rate) of |c), respectively. It can be readily verified that 
the steady solution of Eq. (37) takes the following form 



Pab 



Pcb 



in p [f bc +i{A ab -A ac )] 



2{ (7o6+*A o6 )[ 7i , c +i(A oi ,-A oc )] + 



4 J 



(38) 



l|(7ab+iA ab )[ 7bc + J (A a6 -A ac )] + Si^| 



Apparently, there exists a relation between p c b and p a b, 
i.e., 



Pcb 



9* 



2 {jbc + «(A a6 - A ac ) 



Pab- 



(39) 



Note that since the level pairs |a)-|c) and |a)-|6) can be 
coupled to two laser fields, the parity of level \a) is dif- 
ferent from both \b) and |c). If, for example, \a) possess 
an odd parity, \b) and c) will have an even parity. Thus, 
the electric dipole matrix elements p a b = (a|er|6) =^ 
and p c b — (c|er|6) = 0, and the magnetic dipole ma- 
trix elements m c b = (c|(e/2m e )(L + 2S)|6) ^ and 
m a f, = (a|(e/2m c )(L + 2S)|6) = 0, where L and S denote 
the operators of the orbital angular momentum and spin 
of electrons, respectively. So, it is possible for the nearly 
resonant probe laser to cause the electric-dipole transi- 
tion between levels \a) and \b) as well as the magnetic- 
dipole transition between levels c) and \b) in the three- 
level atomic medium. The electric-dipole transition (|a)- 
|6)) and the magnetic-dipole transition (|c)-|6)) will yield 
the electric polarizability and the magnetic susceptibil- 
ity at probe frequency, respectively. In general, the di- 
mensionless ratio |m c {,/ (p a ;,c) | ~ 10~ 2 in an atomic sys- 
tem, where c denotes the speed of light in a vacuum. 
For this reason, the magnetic-dipole transition may not 
be considered in the treatment for the wave propagation 
in an artificially electromagnetic material. However, in 
a three-level EIT medium where the intensity of cou- 
pling laser is much larger than that of the probe light, 
the population in level |c) is much greater than that in 
the upper level \a). In other words, the stronger cou- 
pling laser enhances the probability amplitude of level 



|c) (according to the relation (39), for the case of small 
detunings, \p cb \ — \(^ c /lbc)Pab\- In a typical EIT ex- 
periment [20,21], njjbc - 10 ~ 100). Thus, thc order 
of magnitude of the density matrix element p c b may be 
larger than that of p a b- This, therefore, means that the 
magnetic dipole moment {2m* cb p c b) may possibly have 
the same order of magnitude as the electric dipole mo- 
ment {2cp* ab pab)- As to the problem of local field cor- 
rection, here, in the simplest case (valid for gases), we 
can take the local field to be the same as the macro- 
scopic field (average field in the sample). So, we need 
not consider the Clausius-Mossotti-Lorentz relation. In 
a three-level atomic system, the electric polarizability 
Xc and the magnetic susceptibility Xm are of the form 
Xc = 2Np* ab p ab / (e £) and Xm = 2Nm* cb p cb /H, respec- 
tively, where N, £ and TL denote the atomic density (to- 
tal number of atoms per volume), the electric and mag- 
netic field envelopes, respectively. Thus, with the help 
of the steady solution (38), one can obtain the relative 
permittivity e r = 1 + Xc and the relative permeability 
/i r = 1 + x m of the above atomic system. By using 
the relation TL = y^e r eo/ p r pa£ between the envelopes of 
magnetic and electric fields, we have the expressions for 
the electric polarizability and the magnetic susceptibility, 



i.e. 



Xc 



Xm 



iN\p ab \ 2 [ lbc +i{A ab -\ aa )} 



toh{ {l ab +iA ab )[ lbc +i{A ab -A a 



m *ab / 1+Xn 
P*a b C V 1 + X<: 



2 7ilc +i(A i,-A ac ) 



(40) 



It follows that the larger is thc Rabi frequency £l c , thc 
more significant is the magnetic-dipole transition caused 
by the quantum interference. So, the magnetic-dipole 
transition deserves consideration in treating the wave 
propagation of the probe light in the case of stronger 
coupling field. 

It is shown that such an EIT medium may have the 
negative permittivity and permeability, and will there- 
fore become an ideal candidate for realizing isotropic left- 
handed media, the optical refractive index of which for 
the probe light is a negative number. Such peculiar me- 
dia exhibit a number of novel electromagnetic and optical 
properties, including reversals of both the Doppler shift 
and Chcrcnkov radiation, anomalous refraction and am- 
plification of evanescent wave [23-26]. Since this subject 
(the scheme of EIT-based realization of negative refrac- 
tive index materials) is beyond the scope of the present 
paper, here we will not consider it further, but will pub- 
lish it elsewhere. 



VI. CONCLUDING REMARKS 

We treat the wave propagation of the probe light in a 
three-level EIT medium and then obtain the solutions of 
Maxwcll-Bloch equations (containing the decay terms) by 
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using the adiabatic approximation under the slowly vary- 
ing envelope condition. The expressions for the time evo- 
lution of both the probe light and the off-diagonal density 
matrix elements are presented, which can be used to con- 
sider the variation of the propagation vector of the probe 
light. This property, i.e., the change of the propagation 
direction and polarization of light in an EIT medium by 
using coherent control may be important for the control- 
lable manipulation of optical fields (e.g., light storage), 
as well as for the quantum information processing. Under 
the conditions of parameters which are taken as in the 
ordinary EIT experiments, the results (e.g., the group ve- 
locity and the electric susceptibility for the probe light) 
obtained here can be reduced to the ones derived from 
the steady solutions of Bloch equations [22] . 
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Appendix 

The three-level density matrix satisfies the following 
equation 



d . 

di P = - 1 
with the Hamiltonian 



H,p 



H 



In —In 



ba 
2^ ca 



2 







U) c 



(Al) 



(A2) 



where Q a b, Sl ac denote the Rabi frequencies of two optical 
fields (the probe and coupling lasers) coupled to the ab 
and ac transitions, respectively. T is a diagonal decay 
matrix of the levels, which agrees with (i|r|j) = juSij. 
uj a , lui, and oj c are the level frequencies of \a), \b) and 
|c), respectively. If the density matrix can be rewritten 
as p(t) = p(t) exp [i (uiijt — kij ■ r)] , then it follows from 
Eq. (A2) that 

Paa = Im (ftpPab + ^cPac) ~ laaPaa, 

pab = -iA ab p ab + | [Q c p cb + Op (fibb ~ Paa)] ~ JabPab, 
p ac = -iA ac p ac + § [fl p pbc + (Pec - Paa)] ~ JacPac, 
pbb = Im (flppba) - IbbPbb, 

pbc = -i (A ac - A ab ) p bc + | (f2 p p QC - Q c pba) - JbcPbc, 
p cc = Im (Cl c Pca) - "fccPcc- 

(A3) 

Here dot denotes the derivative with respect to time. 
fi p (= n a b) and £l c (= fl ac ) are the Rabi frequencies of 
the probe and coupling light, respectively. The frequency 
detunings are defined as follows: A ab = u a — u)b — u a b, 
A ac = Lu a - lu c - uj ac and A bc = A ac - A ab . It 
should be noted that if the probe and coupling light 
fields are in resonance with the ab and ac transitions, 
respectively, namely, the following conditions are satis- 
fied: A a b = A ac = Abc = 0, Eqs. (A3) can be simplified 
into Eq. (1) that has been considered in this paper. 



FIGURE CAPTION 

Fig. 1. The schematic diagram for the A-type EIT sys- 
tem. The level pairs |a)-|6) and |a)-|c) couple the probe 
and coupling fields, respectively. 
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